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Arenaria katoana (Caryophyllaceae) is a rare endemic plant that occurs disjunctly on ultramafic rock 
sites in Japan. This study assessed the degree of genetic distinctness of small, disjunct, isolated popula¬ 
tions of A. katoana in Japan. We sequenced several non-coding regions of the chloroplast DNA ( trnL 
intron, trnL-trnF intergenic spacer, psbA-trnJI intergenic spacer, and petG-trnP intergenic spacer) and 
nuclear ribosomal DNA (ITS1, 5.8S nrDNA, and ITS2) to determine the relationships between popula¬ 
tions and to evaluate the level of intraspecific genetic variation. Eight different chloroplast DNA haplo- 
types were detected, and each population exhibited a unique cpDNA haplotype. Four different ITS types, 
which were shared between populations, were detected. The phylogeny and distribution patterns of the 
cpDNA haplotypes and the ITS types indicated that A. katoana can be divided into three groups; a cen¬ 
tral Hokkaido group (Mt. Yubari, Mt. Bozu, and Mt. Furano-nishi-dake), a Mt. Apoi and Mt. Shibutsu 
group, and a Mt. Hayachine group. The populations of A katoana have been subjected to the effects of 
genetic drift and restricted gene flow, both of which promote population differentiation. 
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Ultramafic areas may have adverse edaphic 
effects on vegetation due to the presence of large 
amounts of Nickel (Ni) and Magnesium (Mg), the 
xeric nature of the environment, and stochastic 
disturbances caused by landslides (Proctor & 
Woodell 1975, Chiarucci et al. 2001, Mengoni et 
al. 2003b, Berglund et al. 2004). At the same 
time, the large number of endemic species that 
occur on such sites suggests ultramafic soils are 
hotspots of biological diversity (i.e., centers of 
floristic differentiation and speciation; Stevanovic 
et al. 2003). 

Arenaria katoana Makino (Caryophyllaceae) 
is a remarkable serpentine alpine plant that is dis¬ 
tributed disjunctly on ultramafic (serpentine and 
peridotite) sites in Japan. It is considered to be a 
serpentine relict species (ultrabasicosaxicolous 
relic; Makino 1908, Kitamura 1950, Toyokuni 
1955). Arenaria katoana var. lanceolata on Mt. 


Apoi, which is known for its large number of en¬ 
demic species, is recognized as distinct from typ¬ 
ical A. katoana on the basis of its narrow leaves 
(Tatewaki 1928). 

Populations of Arenaria katoana, one of the 
rarest species in Japan, have decreased by about 
30% over the last 10 years (Environmental Agen¬ 
cy of Japan 2000) and is thus officially recognized 
as a vulnerable species (VU rank). Arenaria. ka¬ 
toana is currently limited to open, ultramafic 
rock sites located around Mt. Yubari, Mt. Totta- 
betsu, Mt. Apoi, Mt. Furano-nishi-dake, and Mt. 
Bozu in Hokkaido, Mt. Hayachine in Iwate Pre¬ 
fecture, and Mt. Shibutsu and Mt. Tanigawa in 
Gunma Prefecture. The number of individuals in 
each population is very small. 

It is necessary to evaluate the genetic diversi¬ 
ty of endemic serpentine plants with narrow dis¬ 
tributions, as small isolated populations are high- 
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ly susceptible to loss of genetic variation due to 
genetic drift. Information on the level of intra¬ 
specific genetic variation in rare species is vital to 
the identification of appropriate units for conser¬ 
vation (Petit et al. 1998, Newton et al. 1999). 
Moreover, ecological and genetic information on 
serpentine plants can be used to assess how the 
plants have adapted to, and evolved in, such a 
stressful environment. 

In this study, we investigated the genetic vari¬ 
ation and phylogeography of A. katoana using se¬ 
quence information collected from chloroplast 
(cp) DNA and internal transcribed spacer (ITS) 
regions of nuclear ribosomal DNA (nrDNA). Four 
non-coding regions ( trnL intron, trnL—trnF inter- 
genic spacer [IGS], psbA-trnH IGS, and petG - 
trnP IGS) in cpDNA and the internal transcribed 
spacer (ITS) regions ITS1, 5.8S, and ITS2 in the 
nrDNA were used to determine the extent of con¬ 
temporary genetic variation. 

Materials and Methods 

Plant materials 

Forty one leaf samples were collected from 
the six the major remaining localities o fArenaria 
katoana var. katoana on Honshu and Hokkaido 
and A. katoana var. lanceolata on Mt. Apoi, Ja¬ 
pan, in 2004 and 2005. The location of the sites is 
indicated in Table 1 and Fig. 1. 

DNA extraction, PCR amplification and sequenc¬ 
ing 


Total genomic DNA was extracted from sili¬ 
ca-gel-dried leaves (20-50 mg) using a DNA ex¬ 
traction kit (Plant Genomic DNA Extraction 
Miniprep System; Viogene, Sunnyvale, CA) ac¬ 
cording to the manufacturer’s protocol. 

PCR was performed in 20 pL volumes con¬ 
taining 50-100 ng template DNA, 2 pL of 10X 
PCR buffer (TOYOBO, Japan), 0.2 mM of each 
deoxyribonucleotide, 0.2 pM of each of the prim¬ 
er pairs, and 0.5 U of BlendTaq DNA polymerase 
(TOYOBO, Japan). Amplification was performed 
using the GeneAmp PCR System 9700 (Applied 
Biosystems, USA). 

The four non-coding regions the cpDNA, 
were amplified by using the following universal 
primers: (1) “c” and “d” in trnL intron and (2) “e” 
and “f” in trnL—trnF (Taberlet et al. 1991), (3) 
“psbA” and “trnH (GUG)” in psbA-trnH (Hamil¬ 
ton 1999), and (4) “petG” and “trnP” in petG- 
trnP (Hwang et al. 2000). Each PCR cycle fol¬ 
lowed the same program: 30 cycle reaction with 
denaturation at 95°C for 30 s, annealing at 58°C 
for 30 s, and extension at 72°C for 1 min, in addi¬ 
tion to initial denaturation at 95°C for 3 min and 
a final extension at 72°C for 2 min. 

For amplification of the ITS regions, the PCR 
cycle proceeded according to the following pro¬ 
gram: 30 cycle reaction with denaturation at 94° C 
for 1 min, annealing at 56°C for 1 min, and exten¬ 
sion at 72°C for 1 min, in addition to initial dena¬ 
turation at 94°C for 3 min and final extension at 
72°C for 2 min. The ITS region was amplified us¬ 
ing universal primers combination (ITS5 and 


Table 1. Location of study sites including cpDNA haplotypes and ITS types found at each site. Numbers in parentheses indi¬ 
cate number of cpDNA haplotypes in each population (refer to Fig. 1). 


Locality 


Latitude (E) Longitude (N) Altitude (m) 


Number of cp DNA ITS 

individuals analyzed Haplotype type 


Mt. Shibutsu, Gunma Pref. 

36.54 

139.10 

Mt. Hayachine, Iwate Pref. 

39.33 

141.29 

Mt. Apoi, Hokkaido (A. 

42.06 

143.01 

katoana var. lanceolata ) 

Mt. Bozu, Hokkaido 

42.50 

142.12 

Mt. Yubari, Hokkaido 

43.06 

142.14 

Mt. Furano-nishi-dake, 
Hokkaido 

43.18 

142.18 


2228 

8 

A 

a 

1344 

12 

B(8), C(4) 

b 

810 

10 

D(2), E(8) 

a 

790 

3 

G 

d 

1667 

6 

F 

c 

1331 

2 

H 

c 
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(i) Present distribution 



Fig. 1. (i) Distribution of A. katoana. Location names indicate populations analyzed and listed in Table 1. (ii) Distribution of 
cpDNA haplotypes. (iii) Distribution of ITS types. Letters correspond to cpDNA haplotypes and ITS types in Tables 2 
and 3. Numbers above circles and squares correspond to number of individuals analyzed. 


Table 2. Variable sites of aligned sequences of cpDNA regions of trnL intron, trnL—trnF, psbA trnl!, and petG trnP. Num¬ 
bers in parentheses indicate number of repeated mononucleotides. 


Region (length) 

trnL 

(583- 

intron 
-585 bp) 




trnL- 

(323- 

trnF 
-324 bp) 

psbA-trnH 
(313 bp) 

petG-trnP 
(375 bp) 

Variable sites 

2 

2 

3 

4 

4 


1 

1 



3 

8 

2 

2 

2 

6 

6 

0 

5 

cp DNA Haplotype 

4 

8 

0 

2 

4 

9 

4 

1 

0 

A 

A 

T(10) 

A 

A 

C 

C(8) 

T 

G 

A 

B 

A 

T(8) 

A 

A 

A 

C(7) 

T 

T 

C 

C 

A 

T(8) 

A 

A 

A 

C(7) 

G 

T 

c 

D 

A 

T(10) 

C 

A 

A 

C(8) 

T 

G 

A 

E 

A 

T(10) 

A 

A 

A 

C(8) 

T 

G 

A 

F 

T 

T(9) 

A 

A 

A 

C(8) 

T 

T 

A 

G 

A 

T(9) 

A 

A 

A 

C(8) 

T 

T 

A 

H 

T 

T(9) 

A 

C 

A 

C(8) 

T 

T 

A 
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ITS4), and one of the internal sequencing primers 
combinations (ITS5 and ITS2, ITS3 and ITS4) 
for degraded DNA (White et al. 1990). 

The amplified fragments (20 pL) were treated 
with 4 pL of 1/20 diluted ExoSAP (GE Health¬ 
care Bio-Sciences KK) solution to remove excess 
primers and dNTPs. Fragments were sequenced 
using the BigDye® Terminator v3.1 Cycle Se¬ 
quencing Kit (Applied Biosystems, USA) with 
the primers described above. We analyzed the se¬ 
quences using an ABI3100 sequencer with the 
Sequencing Analysis 3.7 program (Applied Bio¬ 
systems, USA). 

Statistical analysis 

Multiple alignments of the sequences were 
obtained using clustal W software (Thompson et 
al. 1994). Each gap of the sequenced regions was 
treated as a fifth character state. A haplotype net¬ 
work was built by using the computer program 
TCS 1.21 and estimation of the 0.95 probability 
limit for parsimony and the network connections 
was obtained (Clement et al. 2000). As the pa¬ 
rameters of genetic diversity and populations dif¬ 
ferentiation, total QiT ) and within-population di¬ 
versity ( hS ), and Gst and Nst, were estimated 
following the methods described by Pons & Petit 
(1995, 1996) using the programs HaploNst (http:// 
www.pierroton.inra.fr/genetics/labo/Software/). 
In the estimated two parameter, G ST is based only 
on haplotype frequencies, while N sr takes into 
consideration genetic similarities between haplo- 
types and establishes the extent to which muta¬ 
tional differences between haplotypes contrib¬ 
utes towards population differentiation. The 
differences in these two parameters were com¬ 


pared by a permutation test using 1000 permuta¬ 
tions (Burban et al. 1999). 

Results 

Sequence variations 

Eight haplotypes (A-H) were detected in the 
cpDNA sequences. The following variations were 
identified: the trnL intron (583-585 bp) contained 
four substitutions and a length polymorphism of 
mononucleotide repeats (poly-T), trn L-trnF IGS 
(323-324 bp) contained one substitution and a 
length polymorphism of mononucleotide repeats 
(poly-C), and psbA-trnH IGS (313 bp) and petG— 
trnP IGS (375 bp) each contained one substitu¬ 
tion (Table 2). In multiple alignments of the trnL 
intron, trnL-trnF IGS, psbA-trnH IGS, and 
petG-trnP IGS, proportions of polymorphic sites 
including gaps were 1.03%, 0.62%, 0.32%, and 
0.27%, respectively. The trnL intron had the high¬ 
est level of genetic variation. 

Four ITS types (a-d) were detected in the ITS 
sequences. The ITS region, which was 623 bp 
long and contained no indels, included ITS1 (247 
bp), 5.8S nrDNA (145 bp), and ITS2 (231 bp). 
ITS1 contained four substitutions, while ITS2 
had two substitutions (Table 3). In multiple align¬ 
ments of ITS1, 5.8S nrDNA, and ITS2, propor¬ 
tions of polymorphic sites were 1.62%, 0.00%, 
and 0.87%, respectively. The ITS region se¬ 
quenced directly from PCR products did not show 
nucleotide polymorphisms within individual 
samples. The 18 sequences of the four cpDNA re¬ 
gions and the ITS were deposited in the DNA 
Data Bank of Japan (DDBJ) under accession 
numbers AB433180-AB433196 and AB433631. 


Table 3. Variable sites of aligned ITS sequences. 


Region (length) 

ITS1 
(247 bp) 



ITS2 
(231 bp) 

Variable sites 


1 

1 

2 

1 

2 


8 

0 

6 

0 

7 

0 

ITS type 

8 

2 

3 

4 

9 

7 

a 

C 

G 

G 

C 

C 

C 

b 

C 

G 

A 

c 

T 

A 

c 

T 

G 

A 

T 

C 

C 

d 

T 

A 

A 

T 

C 

c 
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(i) 



A E D 


(ii) 



a 

Fig. 2. Parsimonious networks of cpDNA haplotypes (i) and ITS types (ii) constructed using TCS program. Small circles in¬ 
dicate missing haplotypes. Solid and open bars represent nucleotide substitutions and length polymorphisms, respective¬ 
ly. 

Geographical distribution and phylogenetic re- Gst and Nst was significant at the 5% level from a 

lationships of the cpDNA haplotypes and the ITS permutation test. 

types Mt. Apoi had the same ITS types (a) as Mt. 

The geographical distribution of the cpDNA Shibutsu, and Mt. Hayachine had one ITS type 
haplotypes and their ITS types are shown in Fig. (b). The populations in central Flokkaido (Mt. 

1. The phylogenetic relationships of the cpDNA Yubari, Mt. Furano-nishi-dake, and Mt. Bozu) 

haplotypes and ITS types are shown in Fig. 2. had phylogenetically close ITS types (c and d). 

Our results revealed that the eight cpDNA haplo- The total diversity {hT) and within-population di¬ 

types (type A-H) were restricted to specific pop- versity ( hS ) of the ITS was 0.897 and 0.089, re¬ 
lations. Mt. Apoi had two cpDNA haplotypes spectively. Gst&n<\Nst of the ITS were 0.897 and 

(D, E), which were phylogenetically close to those 0.906, respectively. The values of the two param- 

found on Mt. Shibutsu (A). Mt. Hayachine also eters were not significantly different {P = 0.458). 
had two haplotypes (B, C) that were phylogeneti¬ 
cally close to each other. The other populations Discussion 

had only one cpDNA haplotype each. The popu¬ 
lations in central Hokkaido (Mt. Yubari, Mt. Fu- Phylogeographical patterns of cpDNA and ITS 
rano-nishi-dake, and Mt. Bozu) had phylogeneti- variations 

cally close haplotypes (F, G, and H). The total Mt. Yubari, Mt. Furano-nishi-dake, and Mt. 
diversity {hT) and the within-population diversity Bozu, which are all located close to the center of 

(hS) of cpDNA was 1.000 and 0.140, respectively. Hokkaido, had phylogenetically close cpDNA 

Gst and Nst of cpDNA were 0.860 and 0.958, re- haplotypes (F, H, and G) and ITS types (c and d). 

spectively. The difference between the values of Mt. Shibutsu and Mt. Apoi had phylogenetically 
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close cpDNA haplotypes (A, D, and E), and the 
same ITS type (a). Mt. Hayachine had phyloge- 
netically different cpDNA haplotypes (B and C) 
and ITS type (b) compared to other regions. Thus, 
A. katoana can be phylogenetically divided into 
three groups: the central Hokkaido group (Mt. 
Yubari, Mt. Bozu, and Mt. Furano-nishi-dake), 
the Mt. Apoi and Mt. Shibutsu group, and the Mt. 
Hayachine group. 

Genetic variation of alpine plants in Japan has 
recently been studied with respect to their indi¬ 
vidual habitats, and the effects of climate change 
on plant distribution and colonization routes have 
been revealed (Fujii et al. 1995, 1996, 1997, Ikeda 
etal. 2006, Ikeda & Setoguchi 2006). These stud¬ 
ies have suggested the importance of the moun¬ 
tains in central Honshu as refugia during the 
Quaternary climatic oscillation, and that alpine 
plants have had their own histories and genetic 
structure. In the case of A. katoana , the geo¬ 
graphical distribution and phylogenetic relation¬ 
ships of the cpDNA haplotypes suggest signifi¬ 
cant phylogeographic structure within the entire 
distribution range (Nst > Gst, P < 0.05). In the 
ITS types, we observed no significant phylogeo¬ 
graphic structure. The populations on Mt. Shibut¬ 
su, which is in central Honshu, had low genetic 
diversity, with only a single cpDNA haplotype. 
The northern populations on Mt. Apoi and Mt. 
Hayachine had two cpDNA haplotypes. The low 
haplotype diversity and the unique cpDNA hap¬ 
lotypes in each population may reflect genetic 
drift due to the effect of the small population siz¬ 
es. The amount of genetic diversity within each 
population is low, and most of the existing diver¬ 
sity can be explained by a between-population di¬ 
versity component. In a previous population ge¬ 
netic and phylogenetic study of the endemic 
serpentine plant Calystegia collina (Convolvula- 
ceae) high genetic differentiation and low gene 
flow between populations was found, not just be¬ 
tween small outcrops but also between popula¬ 
tions occupying the same large outcrop (Wolf et 
al. 2000). Furthermore, populations of another 
endemic serpentine plant, Alyssum bertolonii 
(Brassicaceae), were reported to be established 
by a small number of founders, with each popula¬ 


tion determined to be a distinct genetic entity, re¬ 
flecting a long history of spatial isolation between 
individual populations and serpentine sites (the 
“ecological island model”; Lefebvre & Vernet 
1990, Mengoni et al. 2003a). The current geo¬ 
graphical distribution of A. katoana is extremely 
narrow and restricted to serpentine sites; more¬ 
over, the populations are small in size. Thus, the 
populations of A. katoana have been subjected to 
the effects of genetic drift and restricted gene 
flow, both of which promote population differen¬ 
tiation. 

In terms of conservation genetics, the three 
phylogenetic groups, Central Hokkaido, Mt. Apoi 
and Mt. Shibetsu, and Mt. Hayachine, can be rec¬ 
ognized as independent units for protection. Ad¬ 
ditionally, different populations within each phy¬ 
logenetic group also contain unique cpDNA 
haplotypes. A conservation program for A. kato¬ 
ana should therefore be planned based on this 
phylogeographic information. 
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